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ABSTRACT 

We study a sample of spectroscopic binaries (SBs) in the local solar neighbourhood 
(d ^ 100 pc and My ^ 4) in an attempt to find the distributions of the period, P, the 
primary mass, mi, and the mass ratio q (— m^/mi), as well as the IMF of the local popu- 
lation of field binaries. The sample was collated using available SB data and the Hipparcos 
catalogue, the latter being used for distances and to refer numbers of objects to fractions of 
the local stellar population as a whole. We use the better-determined double-lined SBs (SB2s) 
to calibrate a Monte-Carlo approach to modelling the q distribution of the single-lined SBs 
(SBls) from their mass functions, /(m), and primary masses, mi. The total q distribution 
is then found by adding the observed SB2 distribution to the Monte-Carlo SB1 distribution. 
While a complete sample is not possible, given the data available, we are able to address im- 
portant questions of incompleteness and parameter-specific biases by comparing subsamples 
of SBs with different ranges in parameter space. Our results show a clear peak in the q distri- 
bution of field binaries near unity. This is dominated by the SB2s, but the flat distribution of 
the SBls is inconsistent with their components being chosen independently at random from a 
steep IMF. 

Key words: binaries: general - binaries: spectroscopic - stars: statistics. 



1 INTRODUCTION 

Much has been published on binary statistics and much of this 
has been concerned with the statistics of spectroscopic binaries 
(SBs) (e.g. Goldberg, Mazeh & Latham 2003; Halbwachs et al. 
2003; Boffin, Cerf & Paulus 1993; Boffin, Paulus & Cerf 1992). 
Among the statistics discussed three of the most important have 
been the distributions of period (P), primary mass (mi) and mass- 
ratio (q — 77i2 /mi, where the primary in this paper is always 
taken to be the more massive star). These three parameters in par- 
ticular are important as they present a set of variables that suffice 
to describe the principal properties of a binary system and its evo- 
lutionary path. Many papers have also been published specifically 
just on the mass-ratio distributions of SBs (e.g. Hogeveen 1991; 
Trimble 1974, 1978, 1987, 1990). The question of whether there is 
a peak in the q distribution near q — 1 has been a matter of de- 
bate for some time, on which we hope to shed some light in this 
paper. For instance the recent paper of Goldberg et al. (2003) found 
a distinct bimodal distribution; however, unlike ours, their sample 
was not confined to a particular volume and so might be expected 
to exhibit biases related to this fact. 

The Initial Mass Function (IMF) of binaries is also of in- 
terest and, together with the period distribution, is important for 
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the validation of star formation models. The IMF and P distribu- 
tions are also important for understanding the chemical evolution 
of the Galaxy: interacting binary systems have more complex evo- 
lutionary pathways by which material can be lost to the Interstellar 
Medium (ISM) than single stars, leading, for example, to important 
systematic corrections to predictions of the carbon yield (Tout et al. 
1999). Binary-system population synthesis models may then give a 
more complete idea of the enrichment of chemical elements in the 
ISM. This kind of modelling has been done before for single stars 
but to date studies of this kind for double stars have had to make 
certain assumptions about the P, mi and q distributions and the 
IMF. 

There have also been studies of stars within galactic clusters 
(e.g. Kahler 1999) which are thus magnitude-limited and occupy 
specific limited volumes, as is the case for our sample, but which 
constitute very specific subsets of stars all of a particular age. Clus- 
ter stars are also in the process of diffusing from the cluster at a 
significant rate over time-scales of the order of 10 s yrs. Hence the 
statistics of cluster binaries, while interesting in their own right, are 
not necessarily representative of binaries in general, or field bina- 
ries in particular and we do not consider them further. 

What has not been done before to our knowledge is a study of 
the statistics of a distance and luminosity-limited sample of field bi- 
naries that is, as far as is possible with existing data, evolutionarily 
unbiased and complete. 
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In this paper we derive the distributions for the period P, pri- 
mary mass mi, mass ratio q = m,2/m\ (m,2 being the mass of the 
secondary so that q < 1) and the IMF for a distance and luminosity- 
limited sample of SBs in the local solar neighbourhood d ^ 100 pc 
and Mv ^ 4 and from there are able to make some deductions 
about these same distributions for the general population of field 
binaries. Because our sample is luminosity-limited, the only sys- 
tems we will be missing completely will be those with primary 
mass mi < 1.1 Mq on the Main Sequence. This luminosity cut-off 
makes the sample as unbiased as possible in an evolutionary sense. 

While the sample is certainly not complete to the distance and 
luminosity limits stated, we believe that it is as complete as it is pos- 
sible to make it without seriously compromising the sample size, 
and so without also compromising the conclusions drawn from the 
study. Although the sample is far from complete in absolute terms, 
we are able to make an assessment of the incompletenesses and 
are thus able to attempt to compensate for them. We thus believe 
that our sample is the best approximation to a volume-limited sam- 
ple possible with the data currently available (to be truly volume- 
limited the sample would have to include all objects to the specified 
distance and luminosity limits). 



2 THE SAMPLE 

The spectroscopic binary data for the study was taken from the 
Eighth Catalogue of Orbital Elements of Spectroscopic Binary Sys- 
tems (Batten, Fletcher & MacCarthy 1989), hereafter referred to in 
this paper as the 'Batten' catalogue, supplemented by other data 
of R.F Griffin, both published (see the synopsis paper Griffin 2000 
and others of that series) and unpublished (private communication), 
hereafter referred to as 'RFG' data. The Batten catalogue was until 
recently 1 the most comprehensive catalogue of SBs available, its 
selection criterion simply being to include all SB data available at 
the time of compilation. It thus encompasses a very diverse range 
of motivations for observation from all the many contributors, and 
so will contain a variety of selection effects, which are largely un- 
known apart from the tendency mentioned below to favour shorter- 
period systems. The inclusion of the RFG data was designed firstly 
to increase the size of the available dataset (many more objects 
having been observed since publication of the Batten catalogue in 
1989, many of them by Roger Griffin), and secondly to attempt 
to compensate for the inevitable bias of the Batten catalogue to- 
wards shorter-period systems (inevitable because of the difficulties 
involved in sustaining consistent observing programmes for longer- 
period systems, as Roger Griffin has been able to do). 

The Batten catalogue consists of 1469 SBs and the RFG data 
of 498 SBs. The combined Batten and RFG dataset (after removing 
duplicated objects) was filtered to a distance of 100 pc and limiting 
absolute magnitude of 4 (d ^ 100 pc and Mv ^ 4) by correlating 
entries with the Hipparcos catalogue. The correlation was done by 
Henry Draper catalogue (HD) number, or by coordinates (corrected 
for precession) if no HD number existed, using the Hipparcos par- 
allaxes and apparent magnitudes to calculate distances and abso- 
lute magnitudes. The cutoff for absolute magnitude was chosen in 
such a way that the sample should (to that magnitude) be as com- 
plete and homogeneous as possible, while at the same time having 
a sufficient number of systems to derive reasonable statistics. The 
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Table 1. Period/evolutionary categories. The numbers of SBs are for d ^ 
lOOpcandAfv ^ 4. 



Category 


Period range / days 


SBls 


SB2s 


Total SBs 


(i) 


P > 500 


73 


15 


88 


(h) 


500 > P > 10 


84 


60 


144 


(iii) 


10 > P ^ 1 


62 


63 


125 


(iv) 


P < 1 


7 


7 


14 


TOTALS 




226 


145 


371 



chosen limiting absolute magnitude of 4.0 translates to an appar- 
ent magnitude of 7.5 at 50 pc and to 9.0 at 100 pc. The Hipparcos 
catalogue is complete to mv = 7.3, and in some areas down to 
m v = 9.0 (Perryman et al. 1997; Schroder & Pagel 2003), so all 
the SBs in our initial dataset of 1803 distinct stars that are closer 
than 50 pc (more precisely 46 pc) and brighter than My = 4 should 
have been identified via the comparison with Hipparcos, and a rea- 
sonable fraction of those closer than 100 pc should also have been 
identified. A fainter limit would have created an inconsistent sam- 
ple that, while having data from more systems at close distances, 
would be missing a lot of fainter systems at greater distances; it is 
already clear from the number in our final sample that many of the 
systems in the SB catalogues lie beyond 100 pc. The fainter limit 
would have led to a much steeper fall-off of completeness with dis- 
tance than we find in our chosen sample, creating a sample that 
was harder to analyse. As mentioned above, because our sample is 
luminosity-limited with an absolute magnitude limit of (My ^ 4), 
the only systems we will be missing completely will be those with 
primary mass mi < 1.1 Mq on the Main Sequence. 

When filtered in this way, the sample consists of 371 SBs: 145 
double-lined SBs (SB2s) and 226 single-lined SBs (SB Is). It is this 
sample that we work with in the analysis that follows in the rest of 
the paper. 



3 DISTRIBUTIONS 

3.1 Period (P) distribution 

The period distribution was found directly from the spectroscopic 
binary data. For later use in evolutionary studies the periods were 
divided into the four categories given in TableQ It is intended that 
the four period categories should correspond roughly to the follow- 
ing four evolutionary scenarios: 

(i) P ^ 500 d. Systems that will interact, if at all, at a late 
stage in their evolution, with the primary well on its way up the 
Asymptotic Giant Branch (AGB). 

(ii) 500 d > P ^ 10 d. Systems that interact as the primary 
evolves onto the Red Giant Branch (RGB) or the AGB. 

(iii) 10 d > P ^ 1 d. Systems that interact on the early 
RGB or as the primary leaves the Main-Sequence and enters the 
'Hertzsprung Gap'. 

(iv) P < 1 d. Systems that interact while on the late Main 
Sequence or earlier (contact systems). 

To investigate the behaviour of the observed period distribu- 
tions within different volumes within our 100 pc sample, the frac- 
tions of SBs (of all systems, single and multiple, determined from 
Hipparcos) in the above period categories were determined for vol- 
umes of radius 20 to 100 pc in steps of 2pc for Mv ^ 4. The 
absolute magnitude cutoff (the same as for the rest of the study) 
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avoids the increasing incompleteness that would otherwise result 
with increasing distance (without the cutoff we would include too 
many faint stars, many of which could reasonably be supposed to 
be unrecognised binaries). The results are shown in Fig.Q The fig- 
ure demonstrates well the quality of the data to 30 pc. For volumes 
much smaller than d = 30 pc the data are of little statistical signifi- 
cance due to the low absolute numbers involved (e.g. at 25 pc there 
are 12, 11,8 and systems in P categories (i) to (iv) respectively 
out of a total of 284 Hipparcos objects). The low numbers cause the 
fractions to vary somewhat erratically below 30 pc and so are also 
of little use in determining the trends. The results will be discussed 
in more detail in Section|4| 

Clearly, a significant fraction of the systems are missing, but 
we can estimate the complete fractions for each category by ex- 
trapolation to Opc. Although the data are for SBs, the extrapola- 
tions will necessarily be the fractions of all binary/multiples (with 
respect to all stellar objects, single and multiple) as all binaries will 
be detectable as SBs at Opc. The extrapolation works because the 
further out in the sample we go, the less likely it is that a binary 
system is known as an SB. Therefore, as we go in the opposite di- 
rection, towards Opc, the fraction of systems that are detectable as 
SBs increases, and we get closer to the values of the unbiassed frac- 
tions. Fitting curves to the data from 100 pc to 30 pc, we were thus 
able to determine what the complete fractions would be at d = pc. 

Categories (i) and (ii) were fitted by cubic polynomials and 
category (iii) by a quadratic curve. The choice of curve was deter- 
mined by its shape as well as the closeness of the fit: quadratics for 
categories (i) and (ii) and a cubic for category (iii) all gave curves 
that changed from convex to concave with distance and hence were 
not plausible. A cubic fit may be preferred on the grounds that it 
is what might be expected for number density within a spherical 
volume (the lOOpc limit is not large enough for the thickness of 
the galactic disc to become relevant (Schroder & Pagel 2003), in 
which case a quadratic fit might have been expected instead) but 
they are also more suseptible to errors in extrapolations due to the 
higher order term. The levelling out of the curves at larger distances 
is roughly what one would expect on the basis of the kind of bias 
expected at different distance regimes. At closer distances an ap- 
parent magnitude bias would be expected, as the brighter a sys- 
tem's mv the more likely it would be to be known as an SB, while 
at larger distances one would expect a random selection bias due 
to the effectively random basis on which systems are chosen for 
study, causing the fraction of SBs to tend to a constant value at 
larger distances. For category (iv) it wasn't possible, or realistic, 
to do anything more with the data than to make a linear fit over 
the same distance range, this category being considerably limited 
in what we can do with it by its low absolute numbers. 

The fits were made over as large a range of distances as possi- 
ble before low absolute numbers made the data unreliable, in each 
case from 100 pc to 30 pc. Clearly it is difficult to estimate the er- 
rors on such extrapolations and the fractions are not necessarily go- 
ing to be very accurate. However the fractions given by quadratic 
fits for categories (i) and (ii) instead of the cubic fits used can be 
used to give some idea of the errors involved. In both case the 
quadratic fits had lower constants, D, and hence lower fractions 
at Opc. Taking these differences into account the accuracy of the 
extrapolations is estimated to be of the order of 10%. Details for 
all four categories are given in Table|2| The extrapolations give the 
following estimates for the complete fractions of SBs in the four 
period categories (equal to the value of the constants D in Table|2j: 
0.140, 0.158, 0.0580 and 0.00154, giving a total binary fraction of 
0.36 ±0.11 (for stars brighter than M v = 4). 



Table 2. Polynomial fits for Period categories: Fraction = Ad? + Bd? + 
Cd + D, where d is distance in pc. There are 36 data for each category 
(30 pc to 100 pc in steps of 2pc). The constant D is also equal to the extrap- 
olated fraction at pc. 
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Figure 1. Period distributions of all SBs in the sample for the four period 
categories given in Table [Tl together with extrapolations to Opc. Fractions 
to a particular distance are the fractions of all systems (single and multiple) 
with My ^ 4 in the Hipparcos catalogue to that distance. 

The above completes the more general analysis of the incom- 
pleteness of the sample. Selective incompletenesses will be consid- 
ered in Section|4] 

3.2 Primary-mass (mi) distribution 

Firstly, the Hipparcos parallaxes and apparent visual magnitudes 
of the systems, mv, were used to calculate the absolute visual 
magnitudes, My, The primary masses of the SBs were then esti- 
mated from Mv by correcting for an average contribution of the 
secondary and applying a mass-luminosity relationship. 

The absolute magnitude of the primaries, My 1 , were esti- 
mated using the following magnitude offsets from the absolute 
magnitude of the system, My: 

M Vl = M v 4- 0.50 (for SB2s) (1) 
M Vl = M v + 0.20 (for SBls) (2) 

The offsets are necessary as the presence of even a visually unseen 
companion can make a significant difference to the magnitude of 
the system (see the discussion of £ Aurigae systems below). 

The offset for SB2s was determined by the clear peak of their q 
values (Section f3.3> near 1, with an average q of « 0.84 (see Fig.|5| 
in that section). This means that the average SB2 secondary is about 
40 per cent less luminous than the primary and so would contribute 
about 0.5 mag to the system's My. (It should be noted in passing 
that the individual magnitudes of the SB primaries and secondaries 
are not usually known, even in the case of SB2s. Where the Batten 
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catalogue sometimes quotes two magnitudes these are in fact the 
maximum and minimum magnitudes of the system if the apparent 
magnitude of the system is variable.) 

The offset for SB Is is a mean figure suggested by two consid- 
erations. Firstly, if the luminosity of the secondary were, on aver- 
age, 30 per cent less than that of the primary (equivalent to a magni- 
tude difference between the primary and the system of 0.3) then the 
contrast of the secondary's spectral lines would be enough for them 
to be visible, and the system would be observed as an SB2 rather 
than an SB 1 . Secondly, the depths of eclipse in £ Aurigae systems 
(typically also catalogued as SBs, for example £ Aurigae itself is in 
the Batten catalogue) are equivalent to the secondary contributing 
0.1 - 0.2 mag. to the magnitude of the system. (See, for example, 
the photometry of the January 1989 eclipse of r Persei given in Hall 
et al. 1991 and of the 1988 eclipse of 22 Vulpeculae in Griffin et al. 
1993. It should also be noted that the only reason that £ Aurigae and 
22 Vulpeculae are known as SB2s is because their secondaries are 
visible in the ultraviolet band. For observations of £ Aurigae sys- 
tems in the visible waveband they are usually only seen as SB Is, 
as for example is 22 Vul which is given in the Batten catalogue as 
an SB1.) From these two considerations we therefore adopted an 
offset of 0.2 mag to account for an average SB 1 secondary's con- 
tribution to the luminosity of the system. 

Alternatively, we can argue from the q distribution that results 
from our calculations. From Fig. [5] we see that typically q ~ 0.5, 
with a rather large uncertainty because of the flat distribution. Us- 
ing the same argument as above for SB2s, this would lead to a con- 
tribution of about 0.13 to My, which is (within the uncertainties) 
consistent with the 0.2 offset we have assumed. To test the effect 
of the offset, we have run the calculation below again with a zero 
offset, and the results are qualitatively similar (see Section l3~3l . 

A number of different mass-luminosity relationships were 
used to determine mi from My 1 according to the evolutionary 
status of the primary. The evolutionary status was determined from 
its location on an Hertzsprung-Russell Diagram (HRD) using the 
B—V colour index from the Hipparcos catalogue and the value 
of Mv 1 already calculated. The HRD was divided into a number 
of regions based on the characteristic regions used by Schroder 
(1998) and Schroder & Sedlmayr (2001) (see Fig.|2j. For main 
sequence stars the mass-luminosity relationship used was for 
stars half-way through their H core-burning phase obtained from 
detailed theoretical stellar models computed with the well-tested 
evolution code of Peter Eggleton (Pols et al. 1998). The other 
mass-luminosity relationships are from Schroder & Sedlmayr 
(2001) and Schroder (1998). We thus have the following set of 
equations for mi: 



SB1 Offset = 0.2 



For Main Sequence primaries, both 

B-V < (M Vl + 1.50)/5.16 
& M Vl > -1.50 

and 

B-V < 
& My 1 < -1.50 

use the following mass-luminosity relation 



mi = (3.57 - 1.40M Vl + 0.311M V / 
-0.027M Vl 3 ) M Q 



(3) 
(4) 



(5) 
(6) 



SB2s 



Equation of main— sequence dividing line' 
I M = 5.1550 (B-V) -1.7500 



* Main seq. SB2s 
□ Blue Loop Giant SB2s 
a K Giant clump (1) SB2s 
+ K Giant clump (2) SB2s 
o Lower RGB SB2s 




-1 







B-V 



SB1s 



'Equation of main-sequence dividing line' 
U = 5.1550 (B-V) -1.3000 



x Main seq. SB1s 
□ Blue Loop Giant SB1 
a K Giant clump (1) SB 
+ K Giant clump (2) SB 
o Lower RGB SB1s 




-1 







Figure 2. Evolutionary categories for SB2s and SBls. These are used to de- 
termine which mass-luminosity relationship to use in estimating the masses 
of the primaries. 



For Blue-loop giants: 

B-V > (A/ Vl + 1.50)/5.16 
Mvj < 0.6 
mi = (-0.852M Vl + 2.81) M 

For K giant clump stars (subgroup 1): 

B-V > (M Vl + 1.50)/5.16 
0.6 < M Vl < 0.8 
mi = 1.8 M 

For K giant clump stars (subgroup 2): 

B — V > (M Vl +1.50)/5.16 
0.8 ^ Mv x < 1-0 

mi = (-0.852M Vl + 2.2) M 

For Lower RGB giants: 

B-V > (M Vl + 1.50)/5.16 
M Vl ^ 1.0 
mi = 1.25 M 



(8) 
(9) 
(10) 



(11) 
(12) 
(13) 



(14) 
(15) 
(16) 



(17) 
(18) 
(19) 



(7) 



The resulting distributions of estimated primary masses are 
shown in Fig. [3] 
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Figure 3. mi distributions determined in Section lX2l 

3.3 Mass-ratio (q — 1x12/ mi) distribution 

The q distribution for SB2s was found directly from the observed 
orbital semi-amplitudes, K2 and K\, the period, P, and the eccen- 
tricity, e, via: 

mi, 2 sin 3 i = 1.036 ■ 10 _7 (1 - e 2 ) 3/2 (K 2 + K 1 ) 2 K 2 ,iP (20) 

where mi, 2 have units of solar masses (Mo), K\ i2 nave units of 
kms" 1 and P has units of days (Hilditch 2001, p46). Hence: 

m 2 m 2 sin 3 i K\ .... 

q= = r^- = -jp-. (21) 

mi mi sin i «2 

The q distribution for SB Is however is not so easy to deter- 
mine. The closest one can get to q directly is the following function 
of the mass function, f(m), and the primary mass, mi: 

/("*) q 3 sin 3 i 
mi (1 + q) 2 ' K ' 

The mass function is calculated from the observed period, P, pri- 
mary orbital semi-amplitude, Ki, and eccentricity, e, 

f(m) = 1.036 • 10~ 7 (1 - e 2 ) 3/2 K x 3 P (23) 

using the same units as for equation |20| The primary mass, mi, 
may be determined as in Section l3~2l 

There are a number of different methods for determining the 
underlying q distribution from the /(m)/mi distribution. Two of 
the most commonly used are: 

(1) Richardson-Lucy iterative method (not used in this paper). 
This is a method that was first developed by Richardson (1972) 
for image restoration in optics and then first adapted for astro- 
nomical use by Lucy (1974) for the deconvolution of unknown 
distributions. It has since been used on a number of occasions 



for deconvolving q distributions from observed distributions (e.g. 
Hogeveen 1991). As it is not used in this paper no further details 
of this method shall be given here. 

(2) Monte-Carlo simulation (used in this paper in a refined 
form). This involves calculating the q 3 sin 3 i/( 1 + q) 2 distributions 
from a variety of postulated q distributions and matching them to 
the observed q 3 sin 3 + q) 2 distribution. 

[For other methods see, for instance, Halbwachs (1987).] 

While the Monte-Carlo method might be considered to be 
somewhat unsophisticated, especially given some of the assump- 
tions that have had to be made in the past to make the method work 
(see below), we have been able to introduce a number of constraints 
and tests that we hope improve it somewhat and make it more ro- 
bust. By contrast, methods such as Richardson-Lucy are not as di- 
rect and are also more dependent on initial assumptions. 

3.3.1 A short review of previous Monte-Carlo studies 

This method has been used on many occasions before, but a num- 
ber of restrictions and/or assumptions have always had to be made 
in order to make it work. For example Boffin et al. (1993) restricted 
the sample to 213 spectroscopic binaries with red giant primaries, 
assumed a constant mass for the primaries and also an average 
value of sin 3 i. Trimble's study (Trimble 1990) was also limited 
to a subset of SBs, 164 in this case, most of which were K giants, 
their primary masses being determined from their spectral classes. 

A major problem encountered in Monte-Carlo determinations 
of q distributions is how to take account of the unknown orbital in- 
clination angle i. The effect of any assumed i distribution is made 
more critical by the dependence of /(m)/mi on sin 3 i. Typically 
the unknown i s have been accommodated by assuming an aver- 
age value for sin 3 i. As well as being a bit crude, and so perhaps 
of better use where more rigorous methods may not be possible, 
this method does have other distinct problems. Boffin et al. (1993) 
show that, using this method, it is possible to obtain correct looking 
results for a decreasing q distribution while giving a totally wrong 
result for the case where f(q) oc 1 jq. In fact, for this reason Boffin 
et al. (1993) discard a simple Monte-Carlo approach despite oth- 
ers such as Trimble (1990) and Hogeveen (1991) finding that it 
gives results similar to more sophisticated approaches such as the 
Richardson-Lucy method. Mazeh & Goldberg (1992) also find that 
it produces erroneous results; in their paper they show two graphs 
of simulations with invented q distributions and demonstrate how 
badly they are reconstructed using this method: an even distribution 
is reconstructed as a decreasing distribution with constant gradient, 
and a distribution increasing towards q = 1 is reconstructed as 
an upside-down U-shaped distribution. They then go on to show 
that these results are a consequence of some of the initial assump- 
tions of the method being invalid. Another procedure is to adopt a 
'model-fitting' approach where the probability of detecting a sys- 
tem with a certain inclination, i, is set by theoretical considerations. 
This however usually involves making somewhat ad hoc theoretical 
assumptions. 

Trimble (1974, 1990) took two approaches to the unknown 
sin 3 i values: the direct approach using an average value of sin 3 i, 
and a model-fitting approach that assumed the probability, p, of 
there being an orbit of inclination, i, to be proportional to a certain 
function: p(i) oc sinz in Trimble (1990), as suggested in Halb- 
wachs (1987), and p(i) oc sin 2 i in Trimble (1974). However the 



6 J. Fisher, K. -P. Schroder, Robert Connon Smith 



average sin i values were themselves determined by assuming that 
the probability of detecting a system was proportional to a partic- 
ular function. Halbwachs (1987) used three methods: the average 
sin 3 i method, and two other methods not discussed here: one due 
to Abt & Levy (1985) and the other due to Jaschek & Ferrer (1972). 
None of these methods is entirely satisfactory. 

3.3.2 Refined Monte-Carlo Method ( used in this paper) 

In the present study we attempt to avoid some of the problems pre- 
viously encountered with this method by introducing some of our 
own refinements, ones that are largely made possible by the use of 
the Hipparcos catalogue. For instance, the subset of better-known 
SB2s was used to derive constraints that could then be applied to 
all SBs. For the masses, we were able to go a step further than 
previous studies by using the Hipparcos parallaxes to determine 
absolute magnitudes, and from there obtain the primary masses di- 
rectly from their evolutionary status determined from their position 
on an HRD (as discussed above in section l3~2t . We thus circum- 
vented the dual problems of having to make assumptions about the 
unknown masses and of being restricted to a limited sample of stars 
of a particular luminosity class. The study was hence opened up to 
stars of any evolutionary status, subject only to a limiting absolute 
magnitude of 4 as in the rest of our study. 

We also used a method of random inclinations, i, to avoid hav- 
ing to assume an average value of sin 3 i. The angle i varies between 
0° and 180°, but in practice systems with inclinations near 0° or 
180° will not be easily visible as the components of their radial 
velocity in our line of sight will be too small to measure reliably. 
Systems with i near to 0° or 180° will thus tend to be missed in 
surveys. For our Monte-Carlo procedure we have assumed there- 
fore that i varies between a minimum cutoff angle ao and 90° (suf- 
ficient as i only appears as sin 3 i in the equations). The probability 
of detecting a system is therefore zero for i less than ao and pro- 
portional to sin i between ao and 90°, i.e. p(i) — for i < ao and 
p(i) oc sini for ao ^ i ^ 90°. The reason for the proportionality 
to sin i is that the projection of the SB's orbit onto the line of sight 
is proportional to sin i, and hence so is the probability of observing 
the system as an SB. This means that a initial uniformly random 
variate, x, has to be transformed to i = arccos(l — x) (the mathe- 
matical reasoning behind the transformation is given in section 7.2 
of Press et al. 1993). This is, we think, more reasonable than pre- 
vious assumptions that have been made (and indeed have had to be 
made). 

A problem with the estimated mass procedure as implemented 
is that the masses for the main sequence primaries (that is, the 
greater proportion of them) are taken from half-way through their 
H core-burning phase. In reality the masses will be 'smudged out' 
to either side of these values. To simulate this when calculating the 
/(m)/mi distribution from random inclinations a 'Smudge Fac- 
tor', £, was therefore introduced so that the simulated /(m)/mi 
values are multiplied by a random factor of up to (1 ± Q. 

The values of ao and C, were found by comparing the 
/(m)/mi distribution for SB2s from two different sources: (i) 
from their q values and random i values, ao ^ i ^ 90° for differ- 
ent values of ao and C, (1000 random i values per q value) and (ii) 
from their mass functions, /(m), and estimated primary masses, 
mi. The /(m)/mi distributions were plotted as histograms, 5 bins 
per 0. 1 on the f(m) jn%\ axis (see Fig.[4]for an example of the his- 
tograms from each source for the ao calibration). 

The main effect of the smudging out of the masses on the 
/(m)/mi distribution from (ii) was that the maximum value of 
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Figure 4. Examples of histograms for the ao calibration taken from the two 
sources given in Section l3.3.2l The £ calibration has already been performed 
so that the maximum values of /(m)/mi are the same from both sources. 



/(m)/mi became greater than the maximum theoretical value of 
0.25. (A smaller additional effect was that the slope of the distribu- 
tion towards the maximum f(m) /mi value was slightly shallower 
than that from (i).) From a range of values of ( from 0.0 to 0.5 in 
steps of 0.1, C = 0.20 ± 0.05 was selected as giving the same max- 
imum value of q 3 sin 3 + q 2 ) = f{m)/m\, from (i) as from 
(ii). 

Once C, had been determined in this way the value of ao could 
then be found by comparing the ratios of the first bin (/(m)/mi 
from to 0.025) to the total frequency. By this means the best match 
was found to be for ao = 20.5 ± 1.0°. 

The determination of ao and £ were both cases demonstrating 
how the better known SB2s could be used to calibrate parameters 
for the Monte-Carlo approach to the SB Is. Not only were we able 
to validate methods before using them for the unknown SB 1 q dis- 
tribution but we were able to determine fine-tuning parameters for 
the simulation as a whole. We did then have to assume that the val- 
ues of ao and £ were the same for SB Is as for SB2s. This would be 
reasonable if the SB2s and SB Is have similar selection criteria and 
if the probabilities of detecting a spectroscopic binary as an SB2 or 
an SB 1 are independent. The similar behaviour of the fractions of 
SB2s and SBls in Fig.0in Section|4]gives us some confidence that 
this is indeed the case. 

To perform the Monte-Carlo simulation a variety of plausible 
q distributions were constructed by dividing the range of q values 
from to 1 into ten equal bins and choosing frequencies for each 
bin. The resultant f{m)/m\ distributions were then calculated for 
random values of the inclination, i, from ao to 90°, p(i) oc sini, 
with random values of the smudge factor between 1 — £ and 1 + C- 
For each data point on the q distribution a thousand random val- 
ues of i were used to make the resultant f{m)/m\ distribution 
as smooth as possible. The total frequencies for the constructed 
/(m)/mi distributions rang ed from 2 x 10 4 to 2 x 10 s . The 
/(m)/mi distribution was divided in each case into 25 bins from 
to 0.25 (the maximum value of f{m)/m\ for q ^ 1 being 0.25). 
We then had to assume, albeit from reasonable arguments, the off- 
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Figure 5. Histogram of the q distribution for SB2s and SBls combined, 
SB2s on the bottom (filled-in) and SBls on top (unfilled), for Monte-Carlo 
simulations of stepped SB1 q distributions. The solid line of the SBls is the 
best fit, the broken line is the next best fit. 

set to add to My to determine My 1 for the SBls (see equation [3] 
and the text afterwards). 

The q distributions tried were systematic variations on the fol- 
lowing 'types' : exponential-like functions increasing in frequency 
towards q = 1 but with the rise starting at varying q values, 'hump- 
functions' with a pronounced maximum at varying q values and 
different levels upon either side, and 'step-functions' with the step 
at varying q values and of varying size. As the resultant /(m)/mi 
distributions proved to be rather insensitive to many of the varia- 
tions tried we settled on a step function as being a minimal solution 
requiring the fewest arbitrary assumptions. 

However, given that other authors have found a preference for 
a bimodal distribution of q, with a secondary peak at around q = 
0.2 (e.g. Staniucha 1979), we also present results from peaked q 
distributions to see whether the resulting /(m)/mi distributions 
match observations better or worse than the simple step function. 

The test used to determine the optimal q distribution was to 
compare two ratios: the ratio of the frequency of the first bin (0 < 
/(m)/mt ^ 0.01) to the second (0.01 < /(m)/mt 0.02) and 
the ratio of the frequency of the first bin to the total for the 5th to 
25th bins (0.04 < f(m)/m\ ^ 0.25). The first gives a measure 
of the rise in frequency at the low end of the distribution while the 
second ratio summarises the relationship of the size of the peak 
(invariably near f(m)/m\ = 0) to the rest of the distribution. The 
observed ratios for the SB1 /(m)/mi distribution were 3.250 and 
1.492 respectively. The best fit to these ratios for the q distributions 
tried was for a q distribution consisting of a level portion at a value 
of 20 for < q ^ 0.8 followed by a step down to zero and then 
continuing at zero for 0.8 < q ^ 1. When normalised to the total 
number of observed SB Is this gave frequencies of 28.3, 28.3, 28.3, 
28.3, 28.3, 28.3, 28.3, 28.3, and for the ten equally spaced bins 
from q = to 1 . 

The Monte Carlo SB1 distribution is normalised to the total 
number of SBls observed in our sample, to make it directly compa- 
rable with the observed SB2 distribution. If we then simply add the 
normalised SB1 and observed SB2 distributions with equal weight 
to produce an overall distribution we are assuming that there is no 
obvious bias in our sample towards observing SB2s. Looking at the 
statistics of SB2s and SBls in the Batten catalogue and in the data 
from Griffin, there is if anything a bias towards SBls: our final sam- 
ple contains 61% SBls and 39% SB2s. This lack of bias can prob- 
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Figure 6. Histogram of q distribution for SB2s and SBls combined for 
Monte-Carlo simulations of peaked SB1 q distributions. SB2s are on the 
bottom (filled-in) and SBls on top (unfilled). The solid line of the SBls is 
the best fit, the broken line is the next best fit. 



Table 3. A table of the two different sources used in the paper for the 
histograms of combined SB2 and SB1 q distributions. Given are the (nor- 
malised) frequencies for the 10 bins, q =0 to 1, and a comparison of the 
main features of the distributions. 



Source 


Combined q 


General features 




distribution 




Observed SB2 


29, 29, 30, 32, 31, 


Peak, q = 0.9-1.0, 


+ stepped SB 1 


32, 36, 45, 27, 76 


Plateau, q = 0-0.6, 


Monte-Carlo 




at 0.4 of max. 


Observed SB2 


27, 30, 34, 33, 29, 


Main peak, q = 0.9-1.0, 


+ peaked SB1 


26, 27, 33,41, 89 


Broad peak, q = 0.2-0.4, 


Monte-Carlo 




at 0.38 of max. 



ably be explained by the fact that most of the data requires long- 
term observing programmes that have only been possible where 
telescopes have been dedicated to the programme and the choice 
of which stars to observe has not been limited by the exigencies 
of time allocation committees. Target lists are compiled on the ba- 
sis of detected variability, usually long before it is known whether 
the target is an SB1 or an SB2, and objects are kept on the target 
list until an orbit has been determined. This is certainly true for the 
Griffin data and seems likely to be true for the Batten catalogue as 
well, most of which dates from the era of long-term programmes at 
national or private university telescopes. 

Adding the Monte-Carlo SB1 q distribution to the observed 
SB2 distribution according to this equal weight prescription gives 
the combined SB q distribution shown for the step function in 
Fig.|5] This fi gure also shows the effect of adding the next best- 
fitting SB1 q distribution, the difference being slight as far as 
the overall shape of the total distribution is concerned. The figure 
clearly shows a peak towards 9 = 1. Furthermore, the peak comes 
primarily from the SB2 contribution, derived directly from the ob- 
served data, and so is unaffected by uncertainties in the SB 1 distri- 
bution. Nonetheless, the q distribution is qualitatively similar if a 
zero offset is used in equation [2] instead of 0.2; quantitatively, the 
first four SB 1 bins in the best fit with zero offset have frequencies of 
32.3 instead of 28.3, and all the remaining bins have frequencies of 
16.1. This has the effect, if anything, of accentuating the q — 1 peak 
in the overall distribution. Curiously, the zero offset case (which 
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corresponds to the fainter component of the SB 1 making no signif- 
icant contribution to the total luminosity) puts some stars into the 
q — 0.9 — 1 bin. These systems presumably have evolved primaries 
that are much brighter than their unevolved companions. 

The distribution with zero offset gives a hint of a second peak 
for low q, so it is worth looking at the best-fitting q distributions 
with a peak. These are shown in Fig|5| However, the fit to the ob- 
served /(m)/mi distribution is much less good than for our pre- 
ferred stepped q distribution, the peak for the best fit peaked dis- 
tribution is not very pronounced (the best fit histogram is not very 
different from the one for a stepped distribution in Fig[5}, and the 
peak at q — 1 is still dominant. 

It is also interesting to compare our best-fitting q distribution 
with the one that would be predicted if we took the components 
at random from a steep IMF. This prediction has been made by 
Tout (1991), who considered an IMF steep above 1 Mq but flat for 
smaller masses. For the SB2 distribution, he took the lower mass 
cut-off for both components to be 1 Mq and found a curve that 
rose steeply from q = to q — 1, similar to our SB2 distribution, 
although not so concentrated around q — 1 . For the SB 1 distribu- 
tion he took the same lower mass cut-off for the primary but chose 
0.2 Mq as the cut-off for the secondary. This gave a q distribution 
with a strong peak at q = 0.2 and a curve that dropped smoothly 
to a low value at q = 1 (see Figure 6 of Tout 1991). The joint dis- 
tribution is thus bimodal, similar to the result found by Staniucha 
(1979) illustrated in Figure 1 of Tout (1991). This strongly double- 
peaked distribution is not consistent with our q distribution, which 
is very flat for q < 0.7, even for the peaked distribution that we 
tried (Figs|5|and|6j. We conclude that the components in our sam- 
ple of binaries were not chosen independently and at random from 
the steep IMF that they seem to obey (see discussion in the next 
section). A similar conclusion was reached by Eggleton, Fitchett & 
Tout (1989) for a more restricted sample of visual binaries with two 
bright components. 

FigQdemonstrates that the fractions of SB2s and SB Is of all 
stars (all entries in HIPPARCOS within 100 pc and the same abso- 
lute magnitude limit) behave rather similarly for varying limiting 
absolute magnitudes (the deviation at brighter My for SB2s being 
due to low absolute numbers). This to some degree justifies our 
using parameters derived from SB2s, such as the value of arj, for 
SB Is as well. Note that here, as in the next section, we are dealing 
with My, and not My 1 , as it is the selection biases on the systems 
that we are interested in. 

Given that the SB2s and SB Is seem to have essentially the 
same selection criteria, the two distributions will also be indepen- 
dent of each other, the probability of detecting one set of lines or 
two then depending only upon the detector resolution. In this case, 
detecting one set of lines will be independent of detecting the other 
and it is thus justifiable to add the two q distributions together in 
the way we have done to produce Fig. [5] 



4 DISCUSSION 

So far, we have derived the observed period, primary mass and 
mass-ratio distributions for our distance-limited sample. If there are 
any serious selection effects, however, the true distributions could 
well be significantly different. Selection effects acting on the q dis- 
tribution have already been discussed in detail in Section [33] and 
thus we can be confident of the reality of the peak in the q distribu- 
tion near q = 1. 

There is however another obvious selection effect possibly 
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Figure 7. Fraction of SBs, distance sC 100 pc, for different limiting My, 
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Figure 8. Frequency distribution of masses (in solar masses) of SB pri- 
maries and of single stars. 
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Figure 9. Fraction of SBs, My ^ 4, for different volumes up to 100 pc. 
Fractions are denned as in Fig.lTI 
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Table 4. Comparing the change in fraction with limiting My (Fig.0 with 
the change with distance (Fig.l9*l. 



Fig. |7| change of fraction with Fig. |5| change of fraction with 
My. volume (max. distance of sam- 

pie). 

Decreasing average my & de- Decreasing average ray but 
creasing average mass going same average mass going from 
from left to right. left to right. 

Fraction falls by a factor of ~ 4.3 Fraction falls by a factor of ~ 5 
from M v = 1 to 4. from 25 to 100 pc. 

Could be due to a shallower Could only be due to an increas- 
PDMF and/or a selection effect. ing incompleteness with increas- 
ing vol. & decreasing ray. 

Factor is approximately the same, therefore the decrease with Mv 
(and hence PDMF) is due to a selection effect. Therefore the PDMF 
of SB primaries is the same as that of single field stars. 

acting on the mi distribution that has not yet been considered: the 
possibility of a lower detection rate for less luminous binaries. This 
would be reflected in a less pronounced increase in the observed mi 
distribution towards smaller masses compared to the true present- 
day mass function (PDMF). In the observed distribution in Fig. [8] 
dN/dlogmi oc mj~ 2,8 , while for single stars it is approximately 
tx m -4,8 (Schroder 1998). However, we need to know if this dif- 
ference is genuine or due to a selection effect (or possibly both). To 
do this we look at the variation in detected SB fraction with volume 
(Fig. |3 and compare it with the variation with limiting Mv within 
100pc(Fig.|7j. 

Table [4] summarizes the following argument. The variation 
with Mv shows fractions with decreasing average mass and de- 
creasing average apparent brightness as the absolute magnitude 
limit becomes fainter, while the variation with volume shows frac- 
tions which again have decreasing apparent brightness as the vol- 
ume is increased but now have the same average mass for all vol- 
umes. The variation with volume shows a decrease in fraction by 
a factor of ~ 5 from 25 to 100 pc, while the variation with limit- 
ing Mv shows a decrease by a factor of ~ 4.3 over a corresponding 
range of Mv = 1 to 4 (a factor of 4 in distance being equivalent to a 
difference of 3 in magnitude). The decrease with limiting My could 
again be due to a shallower PDMF or a selection effect, but the de- 
crease with volume could only be due to the increasing incomplete- 
ness as the volume enlarges (and apparent brightness decreases). 
The fact that the two fractions fall off by approximately the same 
factor shows that the decrease with Mv, and hence the shallower 
PDMF, is indeed due to a selection effect. The true PDMF and IMF 
of the binary primaries are therefore nearly identical to those of sin- 
gle field stars in the solar neighbourhood (Schroder & Pagel 2003 
and references contained therein). 
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